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1. Introduction 



Numerical simulations with lattice QCD are bound to rather small lattice volumes when deter- 
mining the hadron spectrum and other low energy parameters in QCD. The computed observables 
show a volume dependence and a thorough understanding of these effects is important for a correct 
interpretation of numerical data. We report on recent progress related to analytical finite volume 
studies in case of the pion mass. 

A long time ago Liischer established an asymptotic formula [jl|] which relates the size depen- 
dence of the pion mass M n with the 7T7r-forward scattering amplitude ^ K (y), 

M *t- M * = - l6K 2 MnL fjy^{iy)e-^y iL + o{e-Vv~^ n L ) 

with M n i the pion mass in finite volume and where the constraint M n L 3> 1 is assumed. The con- 
tributions in the low energy region are enhanced, where the TTTT-scattering amplitude is represented 
very accurately through its chiral representation [||, ||]. At leading order in the chiral counting the 
forward scattering amplitude assumes to be a constant, and the integral simplifies to a modified 
Bessel function of the second kind. When Colangelo and Diirr evaluated subleading finite volume 
effects with Luscher's formula, they observed numerically enhanced corrections with respect to the 
leading order [Qj. These investigations show a clear necessity to go beyond leading order calcula- 
tions in order to have the finite volume effects of the pion mass under control. Furthermore, as the 
derivation of eq. ( |Q| ) takes into account only exponential contributions of the order exp(— M n L) 
and systematically drops those of the order exp (— 2M W L) , the question arises, whether the 
missed terms might turn out to be numerically relevant. The same question also concerns the 
asymptotic formulae of decay constants [Q]. As we will show, a full two-loop calculation of the 
pion mass in finite volume ChPT clarifies these open points. 

The two-loop calculation appears to be interesting also for its own right. To date, a number 
of finite volume calculations have been performed at one-loop order [||], but as far as we know 
two-loop calculations have only been performed for the quark condensate |7|] and for low-energy 
observables in a closely related field [pp. As finite volume effects occur first at the one-loop level, 
only a two-loop calculation - or alternatively if existing, an asymptotic formula a la Liischer - 
leads to a better understanding of the convergence behaviour of the perturbative expansion. 



2. ChPT in finite volume 



Chiral Perturbation Theory (ChPT) is the effective theory for QCD at low energies. It is nowadays 
a mature field which has been applied successfully for a variety of phenomena, in particular in the 
meson sector. For an introduction and a current status of the field, we refer to Ref. [g, |Uj. 

The effective framework is still appropriate, when the system is enclosed by a box of size 
V = L 3 . We refer to the literature for the foundations [ pT| , 12, 13] and a recent review [14]. Here, 
we only remind of the fundamental results which guided the present calculation: the volume has to 
be large enough, such that ChPT can give reliable results, 2F K L 3> 1. The perturbative calculation 
is bound to the value of the parameter M K L. Whether it exhibits to be large (M K L ^> 1, "/^-regime") 
or small [M % L < 1, "£-regime") implies a different power counting. However, in both cases the 
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effective Lagrangian is the same as in the infinite volume. Here, we only cover the "/^-regime", 
where the system is distorted mildly and the only change brought about by the finite volume is a 
modification of the pion propagator due to the periodic boundary conditions of the pion fields 

G{x°,x)= £ G (x°,x + nL), (2.1) 

ael? 

with Gq(x) the propagator in infinite volume. 
3. Pion in finite volume 

The pion mass in finite volume is defined by the pole equation 

G{p L y l =0, for p L = (iM nL ,0) , (3.1) 
where G(p°,p)^ Y is the Fourier transform of the connected correlation function 

dp 



Gt/Apr 1 =M 2 + p 2 -L L (p 2 ), (3.2) 

with Ll(p 2 ) the self-energy in finite volume and M 2 the pion mass in the chiral limit in infinite 
volume. A determination of the pion mass amounts thus to an evaluation of the self-energy in a 



loop expansion. At one-loop order the finite volume corrections have been evaluated in [ 1 1]. Here, 
we discuss the main steps which guided the two-loop calculation. A detailed derivation of the 
results will be given elsewhere []TJ]. It is convenient to write the pion mass in finite volume to two 
loops in the following manner, 

M 2 nL =M\-^-^ , 

Ml = M 2 -I (0) , (3.3) 
where EW, r = 0, 1 , 2, denote the contribution of the self-energy with r propagators in finite volume 



with non- vanishing vector n (cf. eq. (2.1) ). These terms shall be discussed in some detail. 
3.1 Self-energy to O'th order: 

The contributions to are volume independent by definition and merely renormalize the pion 



mass, cf. eq. fl3.3[ ). A detailed discussion of this calculation can be found in \tm, with which we 
agree. 

3.2 Self-energy to l'st order: 

As Liischer showed [[I]], the leading finite volume effects are captured in TP , > which can be summed 
up in closed form 

= 2 J L^^^^^q-p-q) , (3-4) 
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Figure 1: Integration contour in the complex q\ plane (dashed line) with the pole from the pion propagator 
(diamond) and the branch cut from the 7T7r-scattering amplitude. 

with p = (iM K ,0), T n7 i(p,q, —p, —q) the 4-point function of 7T7r-scattering in the forward scattering 
kinematics and m(n) the number of integer vectors z with z 2 = n. Eq. (3.4) may still be simplified 



considerably. We perform a contour integration in the complex plane of the first component q\. 
Singularities are met by the pole of the propagator as well as the branch cuts from the propagator 
and the 4-point function. At the two-loop level only the cuts of the 4-point function have to be 
considered, lying on a hyperbola in the complex q\ plane, as illustrated in fig. [T[ Performing the 
contour integration along the dashed line, we obtain two terms, one from the residuum of the pole 
I p and the other from the integral along the new integration path. The latter contribution vanishes 
as we push the integration lines to infinity, except for the one along the cut, to be denoted by I c in 
the following, 

L {i) =I p + I c . (3.5) 
Along the lines discussed in [Q], the former leads to the resummed asymptotic formula 

oo 



— oo 



with & % {iy) the 7T7r-forward-scattering amplitude. Restricting the sum to the first addend, we 



recover Liischer's formula, eq. (1.1). Its extension to eq. (3.6) has already been suggested in [14, 
[T6[ ]. Numerically, the resummation turns out to be relevant for moderate X n . 

Concerning the contributions of the cut I c , we make use of a dispersive treatment. As these 
terms only start at the two-loop level, they are both suppressed in the chiral as well as in the large 
L expansion and are therefore expected to be small. This is exactly what we observe numerically. 

3.3 Self-energy to 2'nd order: E^ 2 ^ 

Contributions from two pion propagators in finite volume are finally captured in Y,( 2 \ Notice that 
the corresponding Feynman diagrams are uv-finite and need not to be renormalized. Therefore, it 
only remains to find a feasible numerical representation for these terms. 

4. Numerics 

The numerical analysis is performed in line with the setup of Ref. []!(]]. In fig. ||, we evaluate the 
relative finite volume shift 

_ M JzL -M 7t 

Rm « = M n ' (41) 
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Figure 2: Rm % = M % ^/M n — 1 vs. M K for L = 2, 3,4fm. For explanations of the legend, see text. 



for L = 2,3,4 fm as a function of M n . We show the one-loop (LO) as well as the two-loop result 
(NLO). These shall be compared with the resummed asymptotic formula eq. (3i)) with LO/NLO/ 
NNLO input for the 7T7r-scattering amplitude. Note that the one-loop result and the resummed 
asymptotic formula to LO agree with each other. The best estimate for Rm k is finally obtained by 
adding to the asymptotic pure three-loop contribution the two-loop result (NNLO asympt. + NLO 
non-asympt.). At NLO, the finite size effects encounter low energy constants, leading to a non- 
negligible error band which is only shown for the best estimate. We take up a point already alluded 
in the introduction, namely the large contributions when going from LO to NLO in the asymptotic 
formula (dotted to thin-dash-dotted). Compared to this gap, the additional contributions from the 
full two-loop result (thick-dash-dotted) are very small. Consider eg. a pion mass of M n = 250 MeV 
in a 2 fm box. We find Rm k = 0.0236(41), of which 0.0010 stem from the two-loop corrections 
which are not included in the asymptotic formula. The two-loop and the NLO result from the 
asymptotic formula only drift away, when we go beyond the region where the /j-regime can be 
safely applied. 



5. Conclusions 

We have evaluated the finite volume effects of the pion mass to two loops within ChPT in the p- 
regime. The results are compared with a recently proposed extension of the asymptotic formula of 
Liischer. We find that contributions which were neglected in the latter, are numerically very small 
at the two-loop level and conclude that the resummed asymptotic formula is a convenient method 
to evaluate the finite volume effects beyond the leading order. 
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